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Abstract 

A wide variety of complex physical systems described by unitary matrices have 
been shown numerically to satisfy level statistics predicted by Dyson's circular en- 
semble. We argue that the impact of localization in such systems is to provide 
certain restrictions on the eigenvalues. We consider a solvable model which takes 
into account such restrictions qualitatively and find that within the model a gap is 
created in the spectrum, and there is a transition from the universal Wigner distri- 
bution towards a Poisson distribution with increasing localization. 
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A characteristic statistical property of chaotic (as opposed to integrable) states in 
quantum systems is the distribution of their energies. In particular, the nearest-neighbor 
spacing distribution or the long range spectral rigidity of a local set of levels for a wide 
variety of systems in the chaotic regime agree remarkably well with the universal Wigner 
distributions obtained from the Gaussian random matrix theory [|l], The same is also 
true for ergodic quasienergy eigenstates for a variety of periodically driven systems @] 
described by the Flouquet matrix, whose eigenvalues lie on a complex unit circle, and 
belong to Dyson's "circular" ensemble . We will reserve the term Wigner ensemble for 
eigenvalues on the real line. Both ensembles follow the same Wigner distributions in the 
limit of large number of eigenvalues. 

A new problem in this area is the impact of localization on the statistical properties 
of chaotic eigenstates, which leads to deviations from the universal Wigner distributions. 
While attempts have been made to generalize the Wigner ensemble to include such devi- 
ations at a phenomenological level by imposing suitable constraints |^ , it is clear that 
such constraints can not affect the circular ensemble in the same way because the eigen- 
values are already bounded. Nevertheless, numerical studies involving scattering matrix 
for disordered conductors as well as Flouquet matrix for periodically driven systems 
show similar deviations in the spectral properties 0|. It is therefore worthwhile to 
consider an analytic model that can accomodate such deviations in the circular ensemble. 

In this paper, by considering the scattering matrix describing a disordered conductor 
as an example, we will argue that the qualitative effect of localization on the ststistical 
properties of the circular ensemble is to provide certain restrictions on the eigenvalues. 
We will then construct a solvable model that takes into account these restrictions in a 
qualitative way, and show that this leads to a transition in the spectral properties from 
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the universal Wigner distribution towards a Poisson distribution as a function of a single 
parameter related to localization. 

Let us consider a one-dimensional scattering of plane waves of energy E from a po- 



rn is the mass of the incident particle. The 2x2 scattering matrix S has the simple form 



where cos^ = 2ko/ \J Ak^ + k"^ s\v?{ka)\i — kl/K^]^ ip — kga + /i, cos// = cos^ coska. 

The eigenvalues are 6^*^=*='^. In a very crude way, we might think of the case E > Vq 
to mimick a metal, with plane wave states in the region < x < a, while the case 
E < Vq will mimick a finite length insulator with exponentially localized states in the 
region. It is clear that while in the former case the quantity cos 9 can take on all values 
from zero to unity as ko is varied, it becomes restricted to values less than unity in the 
latter case, where k — ip is imaginary and the term k'^ siii^ {ka)[l — k^/k"^] is replaced by 
sinh^(pa)[l + A;Q/p^]. Such a restriction can be interpreted as a constraint on the possible 
maximum oi Tr{S -\- S^) which is proportional to cos^, and the restriction increases 
with increasing "localization" of the waves inside the barrier. In case of a many-channel 
quasi one-dimensional conductor, we can think of the various channels as having different 
incoming energies, and an ensemble of conductors corresponding to different possibilities 
for the values of k. Channels in the metallic regime will correspond to having all possible 
values of 6 and therefore the eigenvalues will be uniformly distributed on the complex 
unit circle without any restriction. On the other hand if the channels are localized, 
the eigenvalues will be distributed in a way consistent with the restriction on the trace 
as mentioned above. This very crude argument suggests that at a phenomenological 
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level, the impact of localization on the eigenvalue distribution of scattering matrices can 
be incorporated by imposing constraints on Tr{S + S*^). This can be done in a way 
suggested by Balian |T0|, namely introducing Lagrange multiplier functions as constraints 



in the joint probability distribution of eigenvalues. In the present work we will choose a 
constraint that has the qualitative features described above, and for which one can, at 
least in principle, solve for all n-point correlation functions of the eigenvalue distribution. 
The hope is that the qualitative effects obtained from such a solvable model will be 
independent of the particular choice of the model. Indeed we will show that the model 
predicts a transition from the highly correlated Wigner towards an uncorrelated Poisson 
distribution in a way that is qualitatively similar to the transition seen numerically for a 
variety of systems. 

For eigenvalues on the complex unit circle, Dyson's circular ensemble is based on the 
basic ansatz of the random matrix theory that for a physical system described by an 
NxN matrix S with eigenvalues e*^",ri = 1,...A^, the joint probabilty distribution for 
the ensemble of all random 5* matrices consistent with given symmetries (unitarity, time 
reversal etc.) can be written quite generally in the form 

m, Om) = n \e''--e''TllHOm). (1) 

m<n m 

Here a is a symmetry parameter and is equal to 1, 2 or 4 for orthogonal, unitary and 
symplectic ensembles respectively. The function w{6) is a Lagrange multiplier function 
which might take care of any system dependent physical constraint [0, and in general 



may depend on various physical parameters. Note that for unbounded eigenvalues of the 
Wigner ensemble, such a constraint is required to keep the distribution normalizable. For 
the circular ensemble the above distribution is already normalizable for w{6) = constant 



and there is in general no need for additional constraint terms. Dyson has shown explicitly 
that the two- level correlation function for the above distribution for w{9) — ^ is identical 
to that of the Wigner ensemble for unbounded eigenvalues in the large N hmit, and 
therefore leads to the same universal Wigner distributions. However, this distribution is 
valid only in the weakly disordered or chaotic regime, and as we argued before, the impact 
of localization can be accomodated phenomenologically by choosing a Lagrange multiplier 
function constraining Tr{S + S^), or equivalently cos^. Because we have no microscopic 
model at this point, we will choose the constraint, with the correct qualitative features, 
such that the model is exactly solvable. 
Our model corresponds to the choice 



Clearly this has the qualitative features mentioned above, where the parameter A will 
serve as a measure of localization; decreasing A increases the constraint on cos 9. We will 
show that this model is solvable in the sense that the spectral correlations can be written 
down in terms of known functions. It turns out that a more general model with two 
independent parameters, which contains our model (2) as a special limiting case, is also 
exactly solvable. Because of its simplicity as well as possible relation to other problems in 
physics, we will start with the more general model, write down the general solution, and 
will come back to our special limiting case when we analyze and interpret the solution. 
The more general two parameter model is defined by the choice 



where we have used the notation {x; g)„ = ]Tk=o{^ ~ ^q'') ■ With the choice a = q^^^, A >> 
1, and q — e~^l^ in the hmit — > oo, or equivalently 5 — > 1, we obtain w{d) = ^2^/'^{l — 
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(2) 




< g < 1, a^g < 1, 



(3) 



COS 6)^^'^ |TT|, which is our model defined in (2). We will first obtain the general solution 
for the model (3), and show that only in the above special limit the impact of localization 
becomes observable in the spectral correlations. In particular we will show that in this 
limit a gap appears in the density. We will also show as an explicit example that in this 
case the number variance obtained from the two-level function shows deviations from the 
Wigner distribution, towards a Poisson limit. Note that in the other limit a = and 
g — > 0, w{9) and the model reduces to Dyson's circular ensemble. 

For simplicity, we will consider only the case where the symmetry parameter a = 2, 
corresponding to the case without time reversal symmetry. We use the method of or- 
thogonal polynomials and write the product term nm<n le*^™ — e*^"| as a Vandermonde 
determinant whose elements form a set of polynomials orthogonal with respect to the mea- 
sure w{6). For our particular choice of w{6) given in eqn. (3), these are the (normalized) 
Szego polynomials generalized by Askey |TT| : 



{q,q)n{q,q) (a^g, q) 

oo 

{a?q,q)n{aq,q)oo{aq,q)c 



1/2 
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_ y {aq;q)k{a,q)n^k{q '^'e^^)^ 
^ {q;q)k{q;q)n-k 

The polynomials satisfy the orthogonality relation 

— / <l>^(e'^ q)^n{e'-^ q)w{d)de = (5) 

ZTT JO 

where the overline denotes complex conjugate. In terms of these polynomials the two-level 
correlation function is given by 

, . ^-1 

K^ie, 0) = ^w{e)^w{(f)) $fc(e*^)$fc(e*<^). (6) 

k=0 
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We now use the unit circle analog of Christoffel-Darboux identity 



J2 ^kizi)^kiz2) 



(7) 



k=o 1 ~ ^2/21 

where Zi = e^^,Z2 = e"^ and we have used the notation = z'^^^niXlz). We obtain 

the large asymptotics of the polynomials by noting that the ratio — 



{m)N-k{aq\q)N 



1 + 0{l/N). Thus for ^ 00, 



fc=0 



a g < 1, 



where in the last line we have used the q-binomial theorem []I2|. The two-level kernel in 



the large N limit can then be written in the general form 



K 



N 



AN-i){e-4,)/2 
2^r 



(g^/^zi, q^/^/z2, aq^''^/zi, aq^''^Z2] q)c 
(gi/V^i, q^^'^Z2, aq^/^zi, aq^l'^/z2] q)c 



1/2 . 



sin[A^(^-0)/2- A] 
[sin(^ - 0)/2] 



(9) 



where the shift A is given by 

A = Im 



In 



{aq^^^Zi,aq^/'^/z2;q)c 



(10) 



{q^/^Zi,q^/^/z2;q)oo 

and we have used the notation {x,y, ..z; q)n = {x; q)n{y', q)n--{z; q)n- For fixed q, in the 
limit 6 ^ (f), this can be simplified and we obtain 



A ^ 2 
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-)Re 



00 „fc+l/2 



A:=0 



Ziq' 



k+1/2 



ae 



-i(9+0)/2 



q 



.k+1/2 



ID 



Writing 1/(1 — xq^^^^"^) = J2'iZo (xq'''^^^'^^ and summing over k first, we obtain the fol 
lowing identity: 



00 qk+1/2 



7^ + 1 ' 



(12) 



,^0. xq-^^i^ l-q—^- ^ i-g' 
The factor (1 — g)/(l — g'+^) — 1 for q « 1, while it is l/(/ + 1) in the limit g — 1. 



In both limits the sum can be explicitly evaluated; it turns out that the result for q 



contains the g << 1 limit, giving a single expression valid for both limits. The result, in 
the limit 6 ^ ({), is 



Eqns. (9) and (13) constitute the solution for large N for the general model defined by 
(3), in the hmit 6 ^ (j). 

We first consider the density of levels given by a{9) = Kn{9,9). Using (6), (9) and 
(13), we get 



Note that the density has a finite N correction to the uniform density N/ 27r of the circular 
ensemble. It is clear that in the N ^ oo limit, the correction might survive only in the 
q ^ 1 limit such that the product (1 — g)A'^ is kept finite. This is precisely the special limit. 



this is indeed the limit where we expect the effect of localization to become observable in 
the spectral correlations. In the rest of our discussions we will restrict ourselves to this 
limit only. 

The expression (14) for the density of levels has one apparently very disturbing fea- 
ture. Although it is properly normalized to A^, the density actually becomes negative for 
sufficiently small values of ^. In fact the condition for the density to remain positive for 
all values of 9 is that the parameter A > Ac — 2N{s/e — l). For 1 << A << Ac, the density 
is positive only for 9 > 9c given by 2\/e^n^sin(6'c/2) ~ 1/A. Thus with decreasing A, 
i.e. increasing localization, 9c increases. We will now show that the negative density for 
A < Ac implies that there exists a gap in the spectrum for 9 < 9c- 

In order to understand the density for A < Ac, we will briefly use an alternative 




(13) 




(14) 



namely q = e 



and a = q^^^ that defines model (2), and as we argued in the beginning. 
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approach based on large "coulomb gas" approximation [0. If we write w{6) = e 
we can interpret the right hand side of eq. (1) as , where the effective 'Hamiltonian' 
H = a J2mj^n sin ^"'^^^ | — J^n^i^n) and the eigenvalues are given by the stationary 
condition 

V'{e) = aP J ^d(j)a{(j)) cot , (15) 

where (t(0) is the density to be evaluated, V is the derivative of V with respect to 9, P 
denotes a principal value integral, and the range / of the integral is determined from the 
normalization Jj d(j)cr{(f)) = N. Expanding cot(y4 — B) and using the normalization, we 
get (« = 2) 

V'i9) = iV cot - + csc^ -P / d<j) , , (16) 

^ ' 2 2 4 '^cotf -cotf ' ^ ' 

where we have allowed for the possibility that the eigenvalues lie in the region |^| > 
dc, Oc < TT. For our model, V{9) = -ln{l - cos6')^/-^ + const.. Using x = [cot f]/[cot ^] 
and y = [cot f ]/[cot y], we can rewrite eq.(16) as 

+ i)t:^ = ^ /" . (17) 

A 1 + te^ J-i X — y 
where we have defined b = cot^(6'c/2), and f{y)dy = a{(f))d(j). This integral can be inverted 



14 to give 



. i^p t , i±i^^. (18) 



The integral can be evaluated explicitly, giving --^=^^^. Going back to the original 
variables, we obtain 



a(^) = ^l±^sin|^cot2|-cot^^, \9\ > 9^. (19) 

The normalization condition gives sin ^ = y+a ~ ^ ^'^^ This agrees with our 

previous result on the existence of the gap as well as its dependence on A. A similar model, 
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with w{6) = has been solved for the density in the saddle point approximation 

in the context of the large N behavior of U {N) lattice gauge theories in two space-time 



dimensions [|1^]. A similar gap was found (at 6 = n), which suggests that the result is 
not peculiar to the particular model we chose; in particular the results from our solvable 
model should be qualitatively valid for models involving qualitatively similar constraints 
on Tr{S + S^). 

The advantage of our solvable model is that we can go beyond the density and eval- 
uate the two-level kernel from which all n-point correlation functions can be calculated. 
However, we can not use (5) and (6) directly because of the gap in the spectrum. The 
existence of the gap suggests that we must allow for this possibility from the beginning, 
and replace eqn. (5) by 

— / <l>Ue''; g)$n(e^^; q)w{e)de = 6m,n (20) 

Although this means that the polynomials are no longer given exactly by (4), we note that 
for small 6c, the density in the large N limit is almost uniform everywhere except near the 
edges. If we restrict ourselves to this uniform density regime, far from the edges, then the 
only real effect of the gap is to affect the normalization. We have taken this into account 
simply by renormalizing the polynomials (4) by a factor ^/C in (20) above. For small 
values of 9c, equivalent to large A, the normalization constant is C ~ ttz:, where c is a 
constant 0(1). We will restrict our following discussions only to the regime 6* ~ tt, where 
the density is approximately uniform, and the kernel KN{6,(j)) becomes translationally 
invariant: 

2iT sm[-2^J 

where we have included the normalization constant C explicitly, and K{(f) — 9) is the 
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complex conjugate of K{6 — (f)). In order to compare with the random matrix theories, 
we have to "unfold" the spectrum by going to a new variable where the mean spacing 
between nearest levels is unity [jl|. This is obtained by choosing the new variables {(, rj) = 
^(1 + j){9, (jy). In terms of these variables the two-level kernel becomes simply 

l^(C-r^)|-C^I^^^^|. (22) 

Note that this looks identical to the two-level kernel of the gaussian random matrix theory 
Q, if we define a new set of variables (C*)''?*) = ^{C^v)- However, in this new variable 
the mean spacing is not unity, but l/C, so the 'unfolding' of the spectrum will take us 
back to the variable {(yV)- 

The two-level kernel can now be used to calculate e.g. the nearest neighbor spacing 
distribution or the long range spectral rigidity. To demonstrate the qualitative effects of 
localization, we will explicitly calculate the number variance for an interval s, defined as 
=< 71^ > — < n >^. Using r = ( — i], this is given in terms of the kernel as ^ 



(SnY = s-2 dr{s ~r)\K{r)\' 
Jo 



s[l -C] + :-[ln{2ns/C) + 7 + 1] + 0(3-'), C = (23) 

where 7 is Euler's constant. As A — > 00, the linear dependence on s cancels exactly and 
the expression reduces to the universal logarithmic dependence on s characteristic of the 
Wigner distribution. However, for any finite A, there is a leftover linear dependence on s 
with the slope increasing with decreasing A (increasing localization). This clearly signals 
a crossover from Wigner towards a Poisson (for which {5nY = s) distribution similar to 
that seen in the case of unbounded eigenvalues [Q, and also similar to the crossover seen in 
numerical studies of the number variance for S-matrix eigenvalues describing transport 

11 



in mesoscopic conductors 0] as well as for A3 statistics (a related measure of the long 
range spectral rigidity of the Flouquet matrix eigenvalues describing time evolution 
of the Fermi-accelerator model 0. Note that if A is related to a physical parameter like 
the conductance which itself scales with A^, then starting from an intermediate case for 
finite N as given in (23), the distribution will scale towards either Wigner or Poisson limit 
depending on whether A scales towards 00 or with increasing N. 

We briefly point out that the general model (3) might include other physically in- 
teresting models. For example in the limit a = and q — > 1^, the function w{6) 
exp[— cos^ |] jl^, which is the model considered in ref. |15 . 

In summary, we constructed a one-parameter solvable model (as a special limit of 
a more general two-parameter solvable model) for the joint probability distribution of 
eigenvalues of unitary matrices which in the large limit leads to a gap in the density. 
The gap increases as a function of the parameter. By analyzing the effect of the gap 
on the number variance, we argued that the model qualitatively describes the effect of 
localization. 
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